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Abstract 

Recently, it has been found that an effective long-range interaction is realized among local 
bistable variables (spins) in systems where the elastic interaction causes ordering of the spins. In 
such systems, generally we expect both long-range and short-range interactions to exist. In the 
short-range Ising model, the correlation length diverges at the critical point. In contrast, in the 
long-range interacting model the spin configuration is always uniform and the correlation length is 
zero. As long as a system has non-zero long-range interactions, it shows criticality in the mean-field 
universality class, and the spin configuration is uniform beyond a certain scale. Here we study the 
crossover from the pure short-range interacting model to the long-range interacting model. We 
investigate the infinite-range model (Husimi-Temperley model) as a prototype of this competition, 
and we study how the critical temperature changes as a function of the strength of the long-range 
interaction. This model can also be interpreted as an approximation for the Ising model on a 
small-world network. We derive a formula for the critical temperature as a function of the strength 
of the long-range interaction. We also propose a scaling form for the spin correlation length at the 
critical point, which is finite as long as the long-range interaction is included, though it diverges in 
the limit of the pure short-range model. These properties are confirmed by extensive Monte Carlo 
simulations. 
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I. INTRODUCTION 



Divergence of the susceptibility is considered one of the characteristics of second-order 
phase transitions. However, it has been pointed out that in spin-crossover type systems 
which belong to the mean-field universality class, the spin configuration is uniform even 
at the critical point. Spin-crossover (SC) materials are molecular crystals, in which the 
molecules can exist in two different states: the high-spin (HS) state and the low-spin (LS) 
state. The HS state is preferable at high temperatures because of its high degeneracy, 
while the LS state is preferable at low temperatures because it has a low enthalpy.- This 
type of competition exists not only in SC materials, but also in charge-transfer materials, 
Prussian-blue type materials, Jahn- Teller systems, and martensitic materials. Such systems 
are generally characterized by the following parameters: the enthalpy difference between the 
HS and LS states, the difference between their degeneracies (or entropies), and the strength 
of the intermolecular interactions. A general classification of types of ordering processes in 
such systems has recently been proposed.- 

An important characteristic of this phase transition is an effective long-range interaction 
caused by an elastic interaction due to the lattice distortion caused by the different sizes of 
the HS (large) and LS (small) molecules, and the spin configuration at the critical point is 
uniform with no large-scale clustering.^ It has also been found that the long-range interaction 
affects dynamical properties.-^»2 In particular, the critical spinodal phenomena predicted by 
the mean-field theory are truly realized. This contrasts sharply with the case of short-range 
models, in which the spinodal phenomena occur as a crossover because nucleation-type 
fluctuations smear out the criticality.— 

This uniform spin configuration is one of the crucial characteristics of the pure elastic 
model without short-range interactions. However, in real materials, we expect that both 
short-range and long-range interactions should exist. For example, if we consider a usual 
Lennard- Jones potential between molecules which depends on the spin states, the model has 
both elastic and short-range interactions.^- In such systems we expect to see ordering clusters 
due to the short-range interaction, though the critical phenomena would still be governed by 
the long-range interaction. Thus, it is an interesting problem to study the crossover between 
short-range and long-range models. In particular, we expect that the correlation length 
of the spin- correlation function is finite in the thermodynamic limit, even at the critical 
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point, as long as any long-range interaction exists. In the present paper, we study how the 
critical correlation length increases and ultimately diverges when the long-range interaction 
vanishes. 

To grasp the general features of the competition between the long-range and short-range 
interactions, in this paper we study a model in which the long-range interactions are those 
of the Husimi-Temperley model, which is the simplest model in which one can study this 
effect. We investigate how the spin-correlation function develops due to the short-range 
interactions. If the long-range interaction is weak, the system shows ordered clusters near 
the critical point of the short-range model, T^^. If we take the cluster size as the unit length, 
the model can be regarded as a pure long-range model, and it shows the critical properties of 
the mean-field universality class at the critical temperature of the model, . This picture 
enables a scaling analysis of the crossover. The difference of the critical temperatures, 
Tc — T^^, is a function of the strength of the long-range interaction. We introduce a formula 
for the critical temperature as a function of the strength of the long-range interaction, and 
we perform Monte Carlo simulations to confirm this relation. A model very similar to ours 
was previously studied by Hastings as a gwell-stirredh approximation for an Ising model on 
a small-world network.— 

A characteristic of the present model is that the correlation length is finite, even at 
the critical point. We study how the cluster size diverges as the strength of the long- 
range interaction decreases, and we propose a scaling form for the divergence, which is also 
confirmed by Monte Carlo simulations. 

II. MODEL : FERROMAGNETIC ISING MODEL WITH NEAREST-NEIGHBOR 
AND WEAK INFINITE-RANGE INTERACTIONS 

A. Hamiltonian 

First, we consider the effects of a weak, infinitely long-range interaction (Husimi- 
Temperley model) on the Ising model with ferromagnetic nearest-neighbor interactions on 
a square lattice. 
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where denotes nearest-neighbor pairs, and cTi = ±1 denotes the Ising spin on lattice 
site i. The critical temperature of this modeU^ is 

2 J 

= — ^ ~ 2.269 ■ ■ • J. (2) 

' ln(l + x/2) ^ ^ 

We adopt the following Hamiltonian for the long-range interaction: 

^HT = -^EE^^^r (3) 

1=1 j=l 

The critical temperature of this modeU^*^ is 

= 4Jo. (4) 

With Jq = J, this critical temperature is equal to that of the mean-field approximation for 
the Ising model on the present lattice. 

For the crossover, we study the following hybrid model, 

H = (1 - a)His + aTiuT, < a < 1 . (5) 

Here, a controls the relative strength of the long-range interaction. 

B. Dependence of the critical temperature on a 

The critical temperature of the model defined by ([5]) changes from Tj^ to T^"^ as a 
changes from to 1. First, we consider the situation in a naive picture. At a temperature 
T, the short-range order is developed by "His, and we assume that Ai'duster spins are tightly 
correlated and behave as one effective spin. In this case, we introduce an effective spin 

{T,},i = 1,...N' = N/N,i^,ter 

-^cluster 

Si= CTj = ^clusterTj, Ti = ±1. (6) 

j'Gcluster i 

Using this effective spin, "Hht is expressed as 

A J ]\t2 N' N' aj Ar N' N' 

J 0^* cluster '±<>'0^ * cluster v;^ v-^ f„\ 

The short-range part has contributions from interactions at the interfaces between clusters, 
and Tiis is given by 

N' 



"His - -Jy A^ciuster ^i'^j ■ (8) 

{id) 



As the clusters grow, the long-range interactions become effectively stronger than the short- 
range interactions, and the critical temperature is given by 

Tc = 4aJoA^cluster- (9) 

If we estimate A^^ciuster using the Ising correlation length C,^^, which has its origin in the 
short-range interaction, it can be written as 

N^lnster ^ ie^y^ = iey-\ (10) 

where rj is the Ising anomalous dimension and the exponent relations are 

a + 2/3 + 7 = 2, (11) 

7 = (2-r7)z/, (12) 

and a = 0, P = 1/8,7 = 7/4? = 1? and 1] = 1/4 in the two-dimensional Ising model.— 
Then, using the relation C,^^ oc (T — T^^)^'^, 

In case a is very small, Tc ~ Tj^, so 

^c-Tf oc(^)^ocai (14) 

This result agrees with Eq. (9) of Hastings' paper.— We can confirm the above picture by 
an exact argument involving the free energy. Let us consider the free energy of the total 
system with a fixed magnetization, m = ^iCTi/N . The partition function is given explicitly 
by 



(15) 

where Z^^(/3,m) is the partition function of the Ising model at the inverse temperature fi 
for a fixed value of m. Therefore, the free energy is given by 

m)/N = \n{Z{/3, m)) = (1 - a)f^{P{l - a), m) - 4a JomV2. (16) 



5 



Here /^^ is the free energy per spin of the Ising model, and we assume that it can be expanded 
around the critical point in the following form 

Thus the critical point of the present model is given by ^-^^Y^|m=o = 0, or 

{I -a) 



AaJo = 0, (18) 



where the susceptibility of the hybrid model ([5]) diverges. If we adopt the relation 



X 



IS 



(T) oc (T - Tf (19) 



the critical point is given by 

Tc - (1 - a)rf oc (4aJo)^(l - a)''^ - (20) 



1 1 



which agrees with flMj) for small a. We note that equation (jT9|) holds only when T is very 
close to TJ^, so equation (120|) holds only when Tc is very close to (1 — a)T^^. Namely, 
equation ( l20l) is only valid for a <C 1. 

For a = 1, 

\imF{/3,m)/N = \im{l- a)f^{f3{l - a),m) -AaJom^ /2 = -m^ - AJotji^ /2 + ■ ■ ■ . (21) 

a— >■! «— j-l 2 

In this case (ITK]) yields Tc = 4Jo, the critical temperature of the Husimi-Temperly model 

To obtain the numerically correct amplitude for Tc{a), we need the Ising susceptibility 
near the critical point for T > T^f^,— 

X^'(/3) = /3Cofiy, 7 = 7/4 (22) 



with Co = 0.962582 • • •, and 



i=^f^ = irt. (23) 



Equation f[T8|) can be written as 

T; / (l-c.)Tf \i / Te-(l-a)Tf - 
Col T, M (l-a)Tjs ^ 



4aJo. (24) 
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We write tc(«) = so 



4 

(l + gv(t^ + a) = (l + tjf_^j . (25) 



Expanding to lowest order in tc and a while setting Jq = J, we get 



tc ^ with A = ( ^^J^ ) ' ^ 1.352745, (26) 



or equivalently, 



2^MF _ 2^18 

This result agrees with ( fT4l) . 



1.773517ai (27) 



C. Monte Carlo study of the a dependence of Tc 

In order to confirm the scaling relation of the previous subsection, we estimated the critical 
temperatures for various values of a by Monte Carlo simulations. Here we fixed both J and 
Jo to 1.0. Therefore, T^^ = 2.269 ■■ ■ and T^'^ = 4 in these units. We used a standard 
Metropolis method, adopting periodic boundary conditions. In most cases, we performed 
500,000 MCS (Monte Carlo steps) for the data with 100,000 MCS for the equilibration. 
From now on, L denotes the linear system size in units of the lattice constant, so the total 
number of spins is L^. 

We estimate a candidate for the critical temperature Tc{a) for each value of a. In order 
to obtain this value, we study the size dependence of the peak position of the so-called 
absolute susceptibility^^ 

X^l^{{M')-{\M\f) (28) 

as a 'critical point' Tc{a,L) for the size L. We expect that the peak position saturates at 
the critical temperature in the thermodynamic limit: 

Tc{a,L) ^Tc{a,(x>). (29) 

In Fig. [H we depict a typical size dependence of the peak for a = 0.001. By a general 
argument we expect the following size dependence: 

Tc(a,L) -Te(a,oo) oc L-^/^ (30) 
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FIG. 1: (Color online) Temperature dependences of x for a = 0.001 and L = 20, 40, 80, 160, and 320 
from below to above. The peak positions are marked by circles. For large systems, they increase 
with increasing L. The left vertical dashed line represents the critical temperature for the pure 
Ising model, and the right-hand line represents the critical temperature for a = 0.001 obtained by 
the Binder cumulant method (see below). 

Here, u is the critical exponent of the correlation length. However, in the present case, the 
critical phenomena belong to the mean-field universality class, and the definition of u is 
subtle. Namely, if we consider the spatial correlation of the Gaussian model, u = 1/2, while 
in the scaling relation in the mean- field universality class, we have the effective v = 2/(i.-'22i'2i 
In the infinite-range (HT) model, distances are not well defined, and only the total number 
of spins, N, has a meaning. Thus, we should rewrite the relation (1501) as 

Tc(a,L)-Tc(a,cx))ociV-^/<^^ (31) 

where we take the latter case [u = 2/d) as we did in a previous paper.- In the present case 
d = 2 and thus N = L"^, which gives 

^^HT = 1, (32) 
which accidentally agrees with that of the short-range Ising model, 

^^is = 1. (33) 

In Fig. [21 we plot the peak position of x by open squares as a function of for several 
values of a. The critical temperature Tc(a,oo) could in principle be estimated by linear 
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FIG. 2: (Color online) Estimates for Tc{a,L) for different values of a vs L ^. a = (pure Ising), 
0.0001, 0.001, 0.01, 0.1, 0.5, and 1 (Husimi-Temperley) from below to above. L = 20, 40, 80, 160, and 
320. The squares denote the peak positions of x, and the circles denote the crossing positions of 
the Binder cumulant for L and L/2. The upper and lower dashed lines mark the exact critical tem- 
peratures for the HT and Ising models, respectively. The inset shows detail for a = 0.001,0.0001, 
and (pure Ising). 

extrapolation in L~^. However, we find a non-monotonic dependence of the peak position as 
a function of for small values of a. (See also Fig. [H) Only when the size becomes large 
enough to show the critical behavior of the HT model, can we apply the scaling relation 
f pOj) . For small sizes, the system behaves like a short-range model, and the peak position 
moves differently. Indeed, we find that in the scaling region, the peak position approaches 
Tc{a, oo) from below. However, for a < 0.01 we find that it decreases with L for small values 
of L. For a = 0.001, we find that the peak position finally increases again when L goes 
from 160 to 320, while for a = 0.0001, it continues to decrease for all values of L considered. 
Thus, we cannot estimate the infinite-system value by a simple extrapolation of the peak 
position in L^^ for a = 0.0001. 

To obtain more accurate estimates for small a, we also estimated Tc{a, L) from the 
crossing point of the fourth-order Binder cumulant,— 

^4(«,L) = 1-^, (34) 
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for L and L/2. When L and a are small, the crossing value of U/^{a,L) is near the Ising 
fixed-point value, Uf-^ ^ 0.61...,^'^ while for larger L and/or a, the crossing moves down 
toward the exact value for the HT model, Uf-'^ = 1 - r^(l/4)/247r2 = 0.27.. .,24-25 where 
r(a;) is the Gamma function. The values of T at the crossing points are shown as circles 
vs for different values of a in Fig. [2J The temperature dependences of f/4(a, L) for 
different a and L are shown in Fig. [31 For a = 0.1, we find the crossing points located near 
Uf^'^, indicating that the critical properties belong to the mean-field universality class. For 
a = 0.01 and 0.001, we find that the crossing points move from near Uf^^ toward Uf^'^ as L 
increases. These results indicate that the critical point of the hybrid model belongs to the 
mean-field universality class for all a > 0. In the case of a = 0.0001, the crossing point of the 
Binder cumulants for {L,L/2) = (320,160) is still near Uf-^. Because we assume that the 
critical behavior for nonzero a belongs to the mean-field universality class, we get a series 
of upper bounds on the critical temperature as the temperature at which U^i^a, L) crosses 
. Lower bounds are given by the cumulant-crossing temperatures. Our best estimates 
for Tc are obtained by linearly extrapolating the crossing temperatures to = . In this 
way, we estimated the T^{a = 0.0001) = 2.281 ± 0.005. In Appendix |A] we show in detail 
how we estimated this value. 

The extrapolated values for Tc{a, oo) — Tj^ are shown on a log-log scale in Fig. HI For 
small a, the data points fall on a straight line of slope 4/7. As mentioned above, to obtain 
more accurate estimates of the critical temperatures from x, we would need to perform MC 
with much larger systems. For small a, the results from the Binder cumulants are in good 
agreement with the power law, Tc(a) - oc a ' . Thus, we confirm the scaling relation 

Tc{a, oo) - Tc(0, oo) oc a^/^^^ _ (35) 

III. CLUSTER SIZE AT THE CRITICAL POINT 

In the pure long-range model, all spins interact with each other. Thus, the concept of 
distance has no meaning, and the system does not show any clustering. On the other hand, 
in the short-range model, the ordering process occurs as a development of short-range order, 
and the cluster size, i.e., the correlation length, represents the extent of the ordering. In 
Fig. [Ht^a) and Fig. Mjo), we depict typical spin configurations at the critical temperature 
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FIG. 3: (Color online) Temperature dependence of the Binder cumulant for (a) a = 0.1, (b) 
a = 0.01, (c) a = 0.001, and (d) a = 0.0001. Points are Monte Carlo data, and the solid lines are 
polynomial fits. The upper and lower horizontal lines are the fixed-point values for the Ising model 
and the Husimi-Temperley model respectively, and the left vertical lines in (c) and (d) represent 
the critical temperature of the pure Ising model. The vertical dashed lines with horizontal error 
bars represent the critical temperatures obtained by extrapolation of the crossing temperatures as 
described in Appendix [Al 

Tc(a) of the short-range model and the long-range model, respectively. A clear difference 
between the two cases is evident. 

Here it should be noted that non- divergence of the correlation length does not mean 
non-divergence of the susceptibility. In the mean-field model, the susceptibility diverges as 
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FIG. 4: (Color online) The a dependence of the normalized critical-temperature difference, [Tc{a) — 
rj^]/[rj^"'^— TJ^] in a log-log plot. The circles and triangles (red online) denote critical temperatures 
obtained from the peak position of x for L = 320 and 80 respectively, and the squares (blue online) 
denote critical temperatures obtained from Binder cumulants. The horizontal dashed line (red 
online) represents y = 1, and the oblique dashed line (blue online) represents the numerically exact 
theoretical estimate, 1.7735170^/"^ (I27p . The latter line, which involves no adjustable parameters, 
agrees very well with the cumulant-generated data for small a. 

\T-T^^\-\ This means that the fluctuation of the magnetization M diverges as 



This fluctuation can be observed as the fluctuation of the uniform density of the spin con- 
figuration. In Fig. [6] we depict typical configurations at T^"^ with different M/N = m. We 
note that the spin configurations are uniform, but the ratio of numbers of up and down spins 
fluctuates. This causes large fluctuations in the magnetization M, but not in the cluster 
size. In the long-range model, large numbers of spins change uniformly. 

In the hybrid model ([S]), the criticality belongs to the mean- field universality class. How- 
ever, short-range order also develops. Thus, we expect a finite correlation length at the 
critical point, which increases as a decreases. In Fig. [71 we depict typical configurations at 
the critical temperature for various values of a. We clearly see that the size of the clusters 
increases with decreasing a. 




(36) 



12 



(a) 
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FIG. 5: (Color online) (a) Spin configurations for Ising model at TJ = 2.269J, L = 100 and (b) 
Husimi-Temperley model at T^'^ = 4Jo. 




FIG. 6: (Color online) Spin configurations for the Husimi-Temperley model for (a)(m) ~ 0.3 at 
Tf^, (b)(m) ~ 0.0 at T^^ and (c) (m) ~ -0.3 at T^^ . 

IV. FINITE-SIZE SCALING OF THE CLUSTER SIZE AT THE CRITICAL 
POINT 

A. Scaling function 

In this section, we study the correlation length at the critical point for several values of 
a. From the relation ( 120|) . we expect the following relation between the correlation length 
S^c at the critical temperature and a : 

ec(«) oc (Te - (1 - «)rf )-is ^ ((4«Jo;^o)V7is)-^^S ^ ^--is/tis (37) 
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(d) 



FIG. 7: (Color online) Typical configurations of the hybrid model at the critical temperature Tc{a) 
for (a) a = 0.0001, (b) a = 0.001, (c) a = 0.01, and (d) a = 0.1. 



for a <^ 1. Moreover, for a ^ 1 we may assume the following finite-size scaling relation 
with the linear dimension of the system L 



Ua,L) = Lf {La^)=Lf [la 



(38) 



where f{x) is a scaling function which is proportional to 1/x for large x, and constant for 
small X. 

In the case of the short-range Ising model, we can estimate the divergence of the corre- 
lation length by making use of the susceptibility: 



X 



NksT 



1 rL 1 



kuT Jo r'^-^+'i 



(39) 



However, in the present long-range interaction model, the value of the magnetization fluc- 
tuates uniformly but not spatially. Therefore, we cannot estimate ^ from x- 
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FIG. 8: (Color online) The disconnected spin correlation function c(r) at the critical point Tc(a) = 
2.332 for q = 0.001. For each value of L, the results of seven independent runs of 10^ MCS each 
are shown. 

B. Measurement of ^ at Tc : Direct measurement of the correlation function 



Here we estimate the correlation length from the spin correlation function c(r) = (cxiaj), 
where r is the distance between the sites i and j, by the following definition 

/o"-/' icir)-ciL/2))rdr 



Sl^'\c{T)-c{L/2))dr 



(40) 



This definition gives the correlation length if c(r) decays exponentially to its large-r value, 
and also for general cases it gives an estimate of the correlation length. In Fig. |H1 we depict 
a typical example of c(r). At large distance, c(r) is constant,'^ proportional toV^- The size 
dependence of ^(a,L) is depicted in Fig. |9l where we confirm that the correlation length 
saturates for large L as expected, even for quite small values of a. The estimated ^(a,L) 
are plotted in the finite-size scaling plot Fig. [101 in which we assume ^ at the critical point 
depends on a as (|55]1 . We find that the data collapse onto a scaling function and thereby 
confirm the theoretical scaling relation (!38l) . 
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FIG. 9: (Color online) The size dependence of the correlation length at the critical point Tc{a). 
Circles, squares, up-triangles, and down-triangles represent T = 2.281 for a = 0.0001, T = 2.332 
for a = 0.001, T = 2.477 for a = 0.01 and T = 2.984 for a = 0.1, respectively. 
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FIG. 10: (Color online) The scaling plot of the correlation length at the critical point. Circles, 
squares, up-triangles, and down-triangles represent T = 2.281 for a = 0.0001, T = 2.332 for 
a = 0.001, T = 2.477 for a = 0.01 and T = 2.984 for a = 0.1, respectively. The linear system 
sizes are L = 20, 40, 80, 160, and 320. The dashed line is proportional to y = 1/x. The data are in 
excellent agreement with the scaling relation (I38p . 
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V. SUMMARY 



We found that in systems with both long and short-range interactions, the long-range 
interaction dominates the critical properties, even if it is infinitely weak. At the critical 
temperature, although the susceptibility diverges, the cluster size does not. At the critical 
temperature, the system has a finite correlation length. 

In this paper, we obtained a formula for the change of the critical temperature as a 
function of the strength a of the long-range interaction, and also a scaling form for the spin 
correlation length at the critical point. 

The crossover of the nature of the order as the length scale changes was studied by a 
Monte Carlo method. We investigated the values of the Binder cumulant at its crossing 
points. It moved from the value of the short-range Ising model to that of the mean-field 
universality class, which enabled us to estimate the critical point systematically. The result 
agrees well with our proposed formula. We further note that our model can be considered 
as a gwell-stirredh approximation for the Ising model on a small-world network.— 

We also proposed a scaling relation for the correlation length at the critical point as 
a function of a. At the critical point, the spin correlation function at large distances is 
constant,- proportional toV^ with a short-range component characteristic of the correlation 
length ^. We obtained the value of the correlation length from the simulated spin correlation 
function, thus providing numerical confirmation of our proposed a scaling function. 

We expect that the results found in this paper are applicable also for real system with 
degrees of freedom corresponding to lattice deformation. A study of such a model will be 
published elsewhere.— We hope this kind of phenomena will be found in future experiments. 
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Appendix A: The estimation of Tc{a) from the crossing points of the Binder cumu- 
lant 

In the present models, the system behaves hke a short-range Ising model for small a 
and L. If we study the crossing point of the Binder cumulants, Ui{a,L), for small values 
of L and at small a, the crossing point gives a value close to that of the Ising model, i.e., 
Uf-^ ~ 0.61 . . .. However, as the size increases, the crossing points approach the fixed-point 
value of the mean-field model, Uf^'^ ~ 0.27. . .. For a = 0.0001, the crossing point of the 
two largest sizes simulated, L = 320 and 160, still stays near 0.57, which is far from Uf^"^. 
Thus, we cannot obtain the critical temperature directly. Here we estimate Tc{a = 0.0001) 
in the following way. We obtain the crossing points for systems with L and L' = L/2 as 
depicted in Fig. [11] Continuous lines for the cumulants as functions of T were obtained as 
polynomial fits to densely spaced Monte Carlo data obtained from simulation runs of up 
to 10^ MCS. We assume the following properties: (1) f/4 at the crossing point for large L 
equals Uf^^, (2) [^(a, L) is a monotonic function of the temperature, and (3) the crossing 
temperature increases monotonically with L. From these assumptions we find in Fig. [TT] 
that Tc{a) is above T = 2.277 which is the crossing temperature for L' = 160 and L = 320. 
Because U4{a,L = 320) crosses U^^"^ at T = 2.289, Tc{a) is below T = 2.289. By linear 
extrapolation with respect to 1/L of the crossing values for L' /L = 80/160 and 160/320 (see 
the inset in Fig. [2]), we estimated the critical temperature as 

Tc(a = 0.0001) = 2.281 ±0.005. (Al) 
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